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The rapidly increasing volume of asteroseismic observations on solar-type stars has revealed a need for automated 
analysis tools. The reason for this is not only that individual analyses of single stars are rather time consuming, but more 
importantly that these large volumes of observations open the possibility to do population studies on large samples of 
stars and such population studies demand a consistent analysis. By consistent analysis we understand an analysis that 
can be performed without the need to make any subjective choices on e.g. mode identification and an analysis where the 
uncertainties are calculated in a consistent way. 

Here we present a set of automated asterosesimic analysis tools. The main engine of these set of tools is an algo- 
rithm for modelling the autocovariance spectra of the stellar acoustic spectra allowing us to measure not only the 
frequency of maximum power and the large frequency separation, but also the small frequency separation and potentially 
the mean rotational rate and the inclination. 

The measured large and small frequency separations and the frequency of maximum power are used as input to 
an algorithm that estimates fundamental stellar parameters such as mass, radius, luminosity, effective temperature, surface 
gravity and age based on grid modeling. 

All the tools take into account the window function of the observations which means that they work equally well 
for space-based photometry observations from e.g. the NASA Kepler satellite and ground-based velocity observations 
from e.g. the ESO HARPS spectrograph. 
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1 Introduction 

Asteroseismology is proving itself as a powerful tool to in- 
crease our understanding of stars. This progress is down due 
to both ground-based facilities such as Elodie (Baranne et 
al. 1996), CORALIE (Queloz et al. 2000), HARPS (Mayor 
et al. 2003), UCLES (Tinney et al.2001) and UVES (Dekker 
et al. 2000) and space-based facilities such as WIRE (Buzasi 
et al. 2000), MOST (Walker et al. 2003), COROT (Appour- 
chaux et al. 2008) and Kepler (Chaplin et al. 2010). And we 
do not expect the fun to stop here. For the future we have 
high expectations regarding the ground-based telescope net- 
work SONG (Grundahl et al. 2009) and the ESA satellite 
PLATO (Catala 2008). 

Solar-type stars are, due to their rich and structural os- 
cillations spectra some of the best targets for asteroseis- 
mology. From the characteristic large frequency separation 
between oscillation modes with different radial order we 
can get a direct measurement of the (acoustic) radius of 
the stars (Christensen-Dalsgaard et al. 2007) and by com- 
bining this measurement with the characteristic small fre- 
quency separation between modes with different angular 
degree we can get estimates of the mass and the age of 
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the stars (Christensen-Dalsgaard et al. 2007). The acoustic 
background that we measure in the stars provides us with in- 
teresting information about activity and convection (Karoff 
2008). For the best observations it is possible to measure 
the frequencies of the individual oscillation modes which 
provide us with an inexhaustible source of information on 
stellar structure and evolution - including: information on 
energy generation and transport, rotation and stellar cycles 
(Karoff etal. 2009). 

2 Methodology 

In this paper we present a set of tools that can be used in 
an automated fashion to perform an asteroseismic analysis. 
The tools consist of the following four main parts: 

I Measure the large frequency separation 
II Measure parameters in the acoustic background 

III Model the autocovariance spectrum {Ai/, Siy, i^max) 

IV Measure stellar fundamental parameters 

The first step performs a search for the characteristic fre- 
quency separation and eventually measures the frequency 
of this large frequency separation (Aiy) and the frequency 
range in which the oscillations modes show a significant 
characteristic spacing (t'max)- The second step uses these 
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Fig. 1 Maximum summed power as a function of fre- 
quency for the simulation of KIC 10454133. The spec- 
trum shows a clear peak just above 2000 /iHz reflecting 
Vniax = 2121.82 ± 4.13/i Hz. The red line shows a fit to 
the spectrum based on a model that includes an exponential 
background and a Gaussian peak. 



first values of Ai' and t'max to measure the amplitude of 
the oscillations. This is done using the method described 
by Kjeldsen et al. (2008). In this step we also obtain mea- 
surements of the characteristic time-scale and amplitudes of 
different components in the acoustic background. The third 
step uses the newly developed technique for modeling the 
autocovariance spectrum of the acoustic spectrum by Cam- 
pante et al. (2010) to obtain more reliable estimates of not 
only Aiy and t^max, but also the small frequency separation 
diy and in the best cases the rotation rate and inclination. All 
these parameters are not only measured in a consistent man- 
ner using the technique developed by Campante et al. (2010) 
they also come with consistent uncertainty estimates. In the 
last and fourth step we use Ai^, i^max and Siy plus parame- 
ters from the Kepler Input Catalog (Latham et al. 2005) as 
input to the SEEK routine developed by P.-O. Quirion in or- 
der to obtain a consistent set of stellar parameters for these 
stars. 

In contrast to most other automated asteroseismic pack- 
ages (Hekker et al. 2010; Huber et al. 2009; Mathur et al. 
2009) the tools presented here all use amplitude spectra cal- 
culated using least squares sine wave fitting (Karoff 2008). 
Generally it only takes of the order of minutes to calculate 
the spectra used in the analysis presented here. 

The tools presented are automated in the sense that they 
can run from the beginning to the end without human inter- 
action, but we do demand that all results are verified by a 
human before a detection is claimed. I.e. if the human eye 
cannot see any significant signal in any of the plots pre- 
sented below a detection is not claimed. We have tested this 
approach on artificial data in the AsteroFLAG framework 
with very positive results (Chaplin et al. 2008). 



3 The large frequency separation 

For the Kepler observations we search for a characteristic 
frequency spacing between 5 and 400 /xHz in the frequency 
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Fig. 2 Maximum summed power as a function of fre- 
quency lag for the simulation of KIC 10454133. The spec- 
trum shows a clear peak just above 50 /iHz reflecting Aiy = 
101.38 ± 2.66/i Hz. The red line shows a fit to the spectrum 
based on a model that includes an exponential background 
and a Gaussian peak. 

interval 20 to 7000 fiHz. This is done by calculating a cor- 
relation matrix of summed power spectra: 

l'=l'max+2.5Al/ 

C(Al/, Z/^ax) = £ PH- (1) 

Here we test 5 radial orders. Of course oscillation modes 
over more radial orders will be excited in many stars, but 
in this first step we measure only the characteristic spacing 
between the 5 most prominent of them. We start by look- 
ing for a frequency range that show signs of a characteristic 
spacing (see Fig. 1). The mean frequency, the width and the 
uncertainty of this frequency range are all found by fitting 
a Gaussian to the one-dimensional correlation matrix as a 
function of frequency, using least squares (see Fig. 1). 

Having a measure of i^max we then measure the charac- 
teristic frequency spacing, i.e. half the large frequency sep- 
aration (Av). This is done the same way as for I'niax, i-e. by 
fitting a Gaussian to the one-dimensional correlation matrix 
as a function of the large separation using least squares (see 
Fig. 2). 

4 Parameters in the acoustic background 

The charateristic time-scale and amplitude of the different 
components in the acoustic spectra are measured using a 
Harvey-like model (Harvey 1985): 

S(^) ^ y ^'^^ + y a,e~^ + c, (2) 
I j 
where a is the amplitude of the component and r is the 
characteristic time-scale. In intensity we generally see activ- 
ity, granulation, faculae, p-mode and chromospheric oscil- 
lations in the Sun and other solar-type stars (Karoff, 2008). 
Of course only activity, granulation and faculae will act 
according to the Harvey-like model in the acoustic spec- 
trum, whereas we assume that the p-mode and chromo- 
spheric oscillations show a Gaussian-like signature in the 
acoustic spectrum (Lefebvre et al. 2008). The spectra are 
smoothed by a Gaussian running mean with a width equal- 
ing two times Av before the amplitudes and characteristic 
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Fig. 3 Autocovariance spectrum of the acoustic spectrum 
of simulation of KIC 10454133 overlaid by the model in 
red. The spectrum shows prominent features at a frequency 
lag close to 100 /iHz equalling the large separation and 50 
/LtHz equalling half the value of the large separation. At 50 
/iHz a clear double peak is seen which is a signature of the 
small separation. 

time-scales are measured. The reason for this is twofold: 
firstly it ensures that we can see components such as faculae 
with small amplitudes compared to the noise and secondly 
it allows us to assume that the noise in the spectra is nor- 
mally distributed so we can use least squares to measure the 
amplitudes and characteristic time-scales. 

5 The autocovariance spectrum of the stellar 
acoustic spectrum 

The autocovariance spectrum of the acoustic spectrum is de- 
fined as a function of the frequency lag L: 

^(^) = ]^ E {P.-P)iPu+L-P), (3) 

where Pi, is the value of the power density spectrum taken 
at the frequency i/, P is the mean value of the power den- 
sity spectrum in the region of interest covering N frequency 
bins. 

In order to model the observed autocovariance spectrum 
we have to build a model of the power density spectrum 
in the region of interest and then calculate the autocovari- 
ance spectrum of this model. This model of the power den- 
sity spectrum is basically identical to the model used by 
Fletcher et al. (2006). As free parameters we thus have the 
large and small frequency separations, one mean rotational 
frequency splitting and inclination, e which takes into ac- 
count the sensitivity to the surface layers, and one mean 
signal-to-noise ratio of the oscillation modes. The modeling 
is done using the general concept of Markov chain Monte 
Carlo as discussed by Campante et al (2010) which allows 
us to provide a robust estimate of the uncertainties. We con- 
volve the power density spectra with the observed window 
functions before we calculate the autocovariance spectrum. 
Also, in order to enhance the visibility of the features in 
the autocovariance spectrum we smooth the power density 
spectrum with a Gaussian running mean with a width of 1 
/iHz (this is done both to the observed spectrum and to the 
model). An example of a modeled autocovariance spectrum 
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Fig. 4 Probability density function of the large frequency 
separation Aj/ for the simulation of KIC 10454133. The red 
region shows the one a most likely solution. 
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Fig. 5 Probability density function of the small fre- 
quency separation parameter for the simulation of 
KIC 10454133. The red region shows the one a most likely 
solution. 

of the acoustic spectrum of a simulation of KIC 10454133 
is shown in Fig. 3. The analysis of the autocovariance spec- 
trum of the acoustic spectrum provides us with estimates of 
not only /S.v, 5v and I'max, but also the mean rotational rate 
and inclination (see Campante et al. 2010 for a discussion of 
the precision with which these parameters can be obtained). 
For the method presented here we are mainly interested in 
estimates of /S.v, 5v and J^max as these are the parameters 
we need for calculating the fundamental stellar parameters. 
Figs. 4 & 5 show the density probability functions for /S.v 
and 5v for the simulation of KIC 10454133. 



6 Stellar fundamental parameters 

We use the SEEK routine developed by P.-O. Quirion to 
estimate stellar parameters based on the asteroseismic pa- 
rameters /S.V, 5v and I'max and stellar atmospheric param- 
eters Tcff, log g and [Fe/H] obtained from the Kepler In- 
put Catalog (Latham et al. 2005). SEEK also utilizes a 
Bayesian approach to provide robust estimates of the most 
likely stellar model and the uncertainties on the stellar pa- 
rameters. These parameters are radius, mass, heavy ele- 
ment abundance, helium abundance, mixing length param- 
eter, luminosity, surface gravity, effective temperature and 
age (Quirion, in prep.). Because the seismic and stellar at- 
mospheric parameters have been obtained in a consistent 
manner the stellar parameters that we obtain are also con- 
sistent and possess consistent uncertainties. 
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Fig. 6 Autocovariance spectrum of the acoustic spectrum 
of the ground-based observations of a Cen A overlaid by 
the model in red. The spectrum shows the signature of the 
small frequency separation as also see in the simulation of 
KIC 10454133. 

7 Application to ground-based observations 

As the tools presented here, especially the one model- 
ing the autocovariance spectrum of the acoustic spectrum, 
takes into account the window function they are well suited 
for analyzing ground-based observations which often suf- 
fer from daily gaps. We have therefore collected time-series 
of high-precision radial-velocities of 10 solar-type stars (a 
Cen A & B, /3 Hydri, Procyon, v Indi, r Ceti, ^ Arae, 7 
Paw, 6 Eri, HD 63077 and t Hor) observed with HARPS 
and UVES from the ESO Advanced Data Products Query 
Form. We are in the process of analyzing these stars in order 
to measure in a consistent manner fundamental stellar pa- 
rameters such as mass, radius, luminosity, effective temper- 
ature, surface gravity and age with consistent uncertainties. 
Fig. 6 shows the modeled autocovariance spectrum of the 
power spectrum of one of these stars a Cen A (Bedding et 
al. 2003) and Figs. 7 & 8 show the density probability func- 
tions for Av and Si' for this star. Using these measurements 
together with stellar atmospheric parameters from Cayrel 
de Strobel et al. (2001) we obtain R/Rq = 1.21 ± 0.03, 
M/Mq = 1.05 ± 0.09, logL/Le = 0.14 ± 0.02, log 
5=4.30±0.01, Tcff = 5760 ± 100 K and age=6.95±2A7 
Gyr for this star which is in general agreement with the 
models made by Miglio & Montalban (2005). 
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Fig. 7 Probability density function of the large frequency 
separation Av from the ground-based observations of a Cen 
A. The red region shows the one a most likely solution. 
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Fig. 8 Probability density function of the small frequency 
separation parameter Dq from the ground-based observa- 
tions of a Cen A. The red region shows the one a most 
likely solution. 
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